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DEGREE SEQUENCE OF RANDOM PERMUTATION GRAPHS 


BHASWAR B. BHATTACHARYA AND SUMIT MUKHERJEE 


Abstract. In this paper we study the degree sequence of the permutation graph G nn asso¬ 
ciated with a sequence n n £ S n of random permutations. Joint limiting distributions of the 
degrees are established using results from graph and permutation limit theories. In particu¬ 
lar, for the uniform random permutation, the joint distribution of the degrees of the vertices 
labelled \nri~\, \nr 2 ~\, ■ ■ ■, \nr s ~\ converges (after scaling by n) to independent random variables 
Di, D 2 ,. ■ ■, D s , where D; ~ Unif(r;, 1 — n), for n £ [0,1] and i £ {1,2,..., s}. Moreover, the 
degree of the mid-vertex (the vertex labelled n/2) has a central limit theorem, and the minimum 
degree converges to a Rayleigh distribution after appropriate scalings. Finally, the limiting degree 
distribution of the permutation graph associated with a Mallows random permutation is deter¬ 
mined, and interesting phase transitions are observed. Our results extend to other exponential 
measures on permutations. 


1. Introduction 


Let [n] := {1,2, • • • , n}, and S n denote the set of all permutations of [n]. For any permutation 
7r n E S n associate a permutation graph G nn = (V(G nn ), E(G nn )), where V(G nn ) = [n] and there 
exists an edge (i,j) if and only if (i — j)(n n (i) — ^n(j)) < 0, that is, whenever i,j determines 
an inversion in the permutation 7r n . The permutation graphs associated with 7r n and tt~ x are 
isomorphic, and the adjacency matrix D Wri associated with the permutation graph G^ n is 




otherwise. 


For a permutation 7r n E S n , define £ n (a) to be the line segment with endpoints (a, 0) and 


(7r n (a), 1), for a E [n]. The endpoints of these segments lie on the two parallel lines y = 0 and 
y = 1 , and two segments have a non-empty intersection if and only if they correspond to an 


inversion in the permutation. Thus, the permutation graph G 7Tn is the intersection graph of the 
segments {£ n (a)}” =1 . Moreover, for every two parallel lines, and every finite set of line segments 
with endpoints on both lines, the intersection graph of the segments is a permutation graph in 
the case that the segment endpoints are all distinct. The associated permutation graph can be 
constructed as follows: arbitrarily number the segments on one of the two lines in consecutive 
order, and read off these numbers in the order that the segment endpoints appear on the other line. 

Permutation graphs were introduced by Pnueli et al. [37] and Even et al. [18], and they showed 
that a graph is a permutation graph if and only if the graph and its complement are transitively 
orientable, that is, an assignment of directions to the edges of the graph such that whenever there 
exist directed edges ( x,y ) and (y,z), there must exist an edge ( x,z ). They also gave a polynomial 
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time procedure to find a transitive orientation when it is possible. Testing whether a given graph is 
a permutation graph can be done in linear time [31]. Permutation graphs are perfect graphs and, 
as a result, several NP-complete problems may be solved efficiently for permutation graphs. For 
this reason, permutation graphs have found applications in many applied problems, like channel 
routing, scheduling, and memory allocation [ ]. Permutation graphs also have applications in 

comparative genomics and bioinformatics [2]. 

Permutation statistics of a randomly chosen permutation, like its cycle structure, inversions, 
descents, and fixed points, are widely studied and are of fundamental importance in the analysis of 
algorithms. For connections of random permutation statistics with determinental point processes 
refer to Borodin et al. [ ]. Recently, Acan and Pittel [ ] studied when cr(n, m), a permutation cho¬ 

sen uniformly at random among all permutations of [n] with m inversions, is indecomposable (refer 
to [7, P ] and the references therein for more on indecomposable permutations). The probability 
p{n, m) that <r(n, m) is indecomposable, is same as the probability that the random permutation 
graph G a t n rn ) is connected. Acan and Pittel [ ] showed that p(n, m) has a phase transition from 0 
to 1 at m n := (6/7r 2 )nlogn. They also studied the behavior of G rJ ( nm ^ at the threshold. 



Vertex Label 


FIGURE 1. Degree distribution of labelled vertices of the permutation graph associated 
with a random permutation of length n = 10 5 . 


One of the most important graph statistic is its degree sequence. Large graphs, like social 
network graphs, are often studied through their degree sequence [35, 36]. For this reason, the 
degree sequence of random graphs has been widely studied over the years. Refer to Chatterjee et 
al. [12] and the references therein for results about random graphs with a given degree sequence. 
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Degree distribution of a random interval graph was studied by Scheinerman [ ]. Diaconis et al. 

[ , ] studied the limits of threshold graphs and interval graphs, and initiated the limit theory of 

intersection graphs. They also considered the degree distribution of threshold graphs [ ]. 

In this paper, using results from the emerging literature on graph and permutation limit theories, 
we study the degree sequence of permutation graphs associated with a sequence n n E S n of random 
permutations. Given a permutation graph G 1Xn denote by d n (i) := Yl'j—i Dn n (i, j) the degree of 
the vertex labelled i G jn], The quantity d n (i)/n will be referred to as the degree proportion of the 
vertex i E [n]. Given a sequence ir n E S n of random permutations, the permutation process is the 
stochastic process on (0,1] defined by 1 

) := MM). (i.i) 

n 

It is shown in Theorem 3.1 that convergence in distribution of the permutation process implies 
joint convergence in distribution of permutations {7r n } ne N hr the sense of Hoppen et al. [26] and 
the degree process 

;= MM). (i.2) 

n 

The limiting degree process can be described in terms of the limit of the permutation process. 
As a consequence, we derive the degree distribution of the uniformly random permutation graph, 
that is, the permutation graph associated with a permutation chosen uniformly at random from 
S n (Corollary 3.3). Figure 1 shows the degree proportion of the vertices in the permutation graph 
associated with a uniformly random permutation of length n = 10 5 . It follows from Corollary 3.3 
that the degree proportion of the mid-vertex, that is, the vertex labelled \n/ 2], converges to 1/2 
in probability (This can also be seen from the fan-like structure in Figure 1 around the point 1/2). 
We show that d n (\n/2~\)/n has a CLT around 1/2 after an appropriate rescaling (Theorem 3.4). 
We also show that the scaled minimum degree in a uniformly random permutation graph converges 
to a Rayleigh distribution with parameter ^ (Theorem 3.7). 

Moreover, we give sufficient conditions for verifying the convergence of the permutation process. 
These conditions can be easily verified for many non-uniform (exponential) measures on permu¬ 
tations. These conditions together with the recent work of Starr [40] can be used to explicitly 
determine the limiting distribution of the degree process for a Mallows random permutation, for all 
/3 E M (Theorem 3.8). For each a E (0,1], the limiting density of d n (\na])/n has a interesting phase 
transition depending on the value of f3: there exists a critical value (3 c (a) such that for (3 E [0, f3 c (a)\ 
the limiting density is a continuous function supported on [a, 1 — a]. However, for (3 > /3 c (a) the 
density breaks into two piecewise continuous parts. If [3 = 1/T denotes the inverse temperature, 
then this is the statistical physics phenomenon of replica symmetry breaking in the low temperature 
regime. 

The rest of the paper is organized as follows: Section 2 contains the basics of graph and per¬ 
mutation limit theories and their connections. Section 3 gives the summary of our main results. 


^Throughout the paper, 7r n will be used interchangeably to denote both the permutation and the permutation process 
depending on the context. In particular, for a E [n] 7r n (a) will denote the image of a under the permutation 7r n . On 
the other hand, for t £ [0,1] n n (t) = n n (\nt\)/n will denote the permutation process evaluated at t. Similarly, d n 
will be used to denote both the degree of a vertex and the degree process. 
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The limiting distribution of the degree processes, relating it to the convergence of the permutation 
process, is in Section 4. Section 5 covers uniformly random permutation graphs, including the CLT 
for the mid-vertex. Conditions to verify the convergence of the permutation process are given in 
Section 6. These conditions are also verified for a general class of exponential measures on permu¬ 
tations. The application of these results to the Mallows random permutation is given in Section 
7. Section 8 contains the asymptotics for the minimum degree of a uniformly random permutation 
graph. 


2. Graph and Permutation Limit Theories 


2.1. Graph Limit Theory. The theory of graph limits was developed by Lovasz and coauthors 
[9, 10, 30], and has received phenomenal attention over the last few years. Graph limit theory 
connects various topics such as graph homomorphisms, Szemeredi’s regularity lemma, quasirandom 
graphs, graph testing and extremal graph theory, and has even found applications in statistics and 
related areas. For a detailed exposition of the theory of graph limits refer to Lovasz [30]. 

Here we mention the basic definitions about convergence of graph sequences. If F and G are two 
graphs, then define the homomorphism density of F into G by 


t(F,G) : = 


| hom(F, G)| 

|F(G)|I^)I, 


where |hom(i ? , G)| denotes the number of homomorphisms of F into G. In fact, t(F,G) is the 
probability that a random mapping <f> : V ( F ) —> V ( G ) defines a graph homomorphism. The basic 
definition is that a sequence G n of graphs converges if t(F , G n ) converges for every graph F. 

There is a natural limit object in the form of a function W G tF, where W is the space of all 
measurable functions from [0, l] 2 into [0,1] that satisfy W(x,y) = W(y,x) for all x,y. Conversely, 
every such function arises as the limit of an appropriate graph sequence. This limit object deter¬ 
mines all the limits of subgraph densities: if H is a simple graph with V(H) = {1, 2,..., \V(H)\}, 
let 


t(H, W) = / TT W(xi,xj)dxidx 2 ■ ■ ■dx\ V (H)\- 

A sequence of graphs {G n } n >i is said to converge to W if for every finite simple graph Ft, 


lim t(H,G n ) = t(H,W). (2.1) 

n— 

These limit objects, that is, elements of W, are called graph limits or graphons. A finite simple graph 
G on [n] can also be represented as a graphon in a natural way: Define f G (x , y) = 1 {( [~nx~|, \ny \) G 
E(G)}, that is, partition [0, l] 2 into n 2 squares of side length 1/n, and define f G (x, y ) = 1 in the 
(i,j)-th square if (i, j) G E{G) and 0 otherwise. Observe that t{H,f G ) = t{H,G) for every simple 
graph Ft and therefore the constant sequence G converges to the graph limit f G . 

The notion of convergence in terms of subgraph densities outlined above can be captured by the 
cut-distance defined as: 


da(f,g)-= sup 

S,XC[0,1] 


' SxT 


[f{x,y) ~ g(x,y)\dxdy 
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for /, g E W. Define an equivalence relation on W as follows: / r\_/ g whenever f{x,y)=g a {x,y):= 
g(ax,ay), for some measure preserving bijection a : [0,1] i-A [0,1]. Denote by g the closure of the 
orbit g a in (W, do)- The space {g : g E of closed equivalence classes is denoted by W and is 
associated with the following natural metric: 

fa(/, 5 ) : = inf dn{f,9o) = inf d a (f a ,g) = inf d{f ai ,g U2 ). 

(7 (T 71,(72 

The space is compact [9], and the metric fa is commonly referred to as the cut-metric. 

The main result in graph limit theory is that a sequence of graphs {G n } n >i converges to a limit 
W E W in the sense defined in (2.1) if and only if fa {f Cn , W) -> 0 [9, Theorem 3.8]. More generally, 
a sequence {W n } n >i converges to W E W if and only if fa (W n , W) 0. 


2.2. Permutation Limits. Analogous to the theory of graph limits Hoppen et al. [25, 26] devel¬ 
oped the theory of permutation limits. For ir n E S n and a E S a , o is a sub-permutation of Tr n if there 
exists 1 < h < ■ ■ ■ < i a < n such that such that fax) < fay) if and only if n n (i x ) < TT n (i y ). For 
example, 132 is a sub-permutation of 7126354 induced by i\ = 3, = 4, *3 = 6. Sub-permutations 

are often referred to as patterns and their combinatorial properties are widely studied (refer to 
Bona [ 8 ], the recent paper of Janson et al. [27] and the references therein). We follow the setup of 
Hoppen et al. [26] (see also Krai’ and Pikhurko [ ] and Glebov et al. [22]), and refer to them as 

sub-permutations because of its similarity to the notion of sub-graphs, the counterpart object in 
graph limit theory. 

The density of a permutation a E S a in a permutation 7 r n E S n 


is 


d(a, ir n ) = \ 


-1 


#{<r E S a : <t is sub-permutation of 7r n } 
0 


if 

if 


a < n 
a > n. 


An infinite sequence (7r n ) ne ^ of permutations is convergent as n —> 00 if d(a , 7 r n ) converges for 
every permutation a. 

Every convergent sequence of permutations can be associated with an analytic object, referred to 
as permuton , which is a probability measure v on ([ 0 , l] 2 , ^([ 0 , l] 2 ) with uniform marginals, where 
.^([0, l] 2 ) is the sigma-algebra of the Borel sets of [0, l] 2 . For an integer n, sample n independent 
points (xi, r/i), (X 2 , 2 / 2 ), ■ • •, (x n , 2 /n) in [0, l ] 2 randomly from the measure u. Let a x and a y be the 
permutations of order n such that x^i) < x a ^ 2 ) < ■■■x<j x (n) and y ay{l) < 2 /^( 2 ) < ■■■y< Ty {n), 
respectively. Define ay 1 o a x as the v-random permutation of order n. For a permuton u and a 
permutation a E Sk, define d(a, v) as the probability that a u-random permutation of order k is a. 
The main result in permutation limit theory [26] is that for every convergent sequence (TT n )n£N of 
permutations, there exists a unique permuton v such that 


d(a, u) = lim d(a,n n ), ( 2 - 2 ) 

n—>■ 00 

for every permutation a. The permuton v is defined as the limit of the sequence (vr^fagi^. On the 
other hand, a sequence of z/-random permutations ( 7 T n ) ne pj converges to v in the limit. 

As in graph limit theory, the above notion of permutation convergence can be metrized by 
embedding all permutations to the space probability measures A4 on [0, l ] 2 with uniform marginals, 
equipped with any metric which induces the topology of weak convergence. To this end, define for 
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any 7 r n E S n , a probability measure v^ n E Ad as: 

di^iTn '■= U n {x,y)dxdy , where U n (x,y) = nl{(x,y) : 7 r n ([nxJ) = Lra/J}, (2.3) 

is the density of v^ n with respect to Lebesgue measure. Like in graph limit theory, o 7Tn has the 
following interpretation: Partition [0, l] 2 into n 2 squares of side length 1/n, and define / % (x, y) = n 
for all (x,y) in the (z, j)-th square if tt n (i) = j and 0 otherwise. The measure o- Kn will be referred 
to as the permuton associated with ir n . 

The space Ad equipped with the topology of weak convergence of measures is compact. Moreover, 
the map it —> is 1-1 and so all finite permutations are contained in Ad. Analogous to the result in 
graph limit, Hoppen et al. [26, Lemma 5.3] showed that a sequence of permutations ir n convergences 
in the sense of ( 2 . 2 ) if and only if the corresponding sequences of measures v nn converges weakly 
to v E Ad. More generally, any sequence of measures {v n } hr Ad converges weakly to a measure v 
if and only if d(a, v n ) — > d(o , v) for all permutations a. 


2.3. Limit of Permutation Graphs. Diaconis et al. [15, 16] studied the limits of threshold 
graphs and interval graphs. They also mentioned that their methods can be used for other classes 
of intersection graphs, which include permutation graphs [16]. However, instead of describing the 
limit object as symmetric function from [ 0 , l ] 2 to [ 0 , 1 ], they represented the graph limit as a 
measure on [ 0 , l ] 2 with uniform marginals. 

To this end, for every measure v on [0, l ] 2 with uniform marginals, they defined a unique graph 
limit object W v E W by specifying t(F, W u ) for all graphs F as follows: 

t(F, W u ) := E J] K ( x i> x j) 

(iJ)eE(F) 

= TT K(xi,Xj)du(xi)du(x2)---du(x\ V f F )i), (2.4) 

/[n -|12|V(F)| - L - L V ' 

where X\,X 2 ,..., X n are independent and identically distributed from u and K : [0, l ] 2 x [ 0 , l] 2 ^ 
[0,1] given by K((ai,bi), (a 2 , b 2 )) = l{(ai - a 2 )( 6 i - b 2 ) < 0}. 

In the graph limit literature it is usually convenient to represent a graphon by a functional 
W : [0,1 ] 2 —y [0,1]. However, Diaconis et al. [16] described permutation graph limits in terms 
of a permuton v. In that case every probability measure u on [0, l ] 2 defines a graph limit W v by 
(2.4). Diaconis et al. [ ] showed that every permutation graph limit may be represented in terms 

of a measure v on ([0,1] 2 ,^([0, l] 2 )) via (2.4) with the two marginal distributions of v both being 
uniform on [ 0 , 1 ]. 

Glebov et al. [ ] pointed out that if (n n )n£N is a convergent sequence of permutations (in the 

sense of Hoppen et al. [26] described above), then the sequence of permutation graphs (G'^- II ) ne N is 
also convergent. Therefore, each permuton v can be associated with an equivalence class W u E W. 
However, the map v —> W v is not one-to-one, as can be seen by the following example: Suppose 
v E A4 be a permuton which is not exchangeable, i.e. if ( X , Y) has distribution u then (Y, X ) has a 
distribution y / v. If 7 r n be a sequence converging to u, then converges to y. Since the graphs 
G nn and G n -1 are isomorphic, W u are are identical, but by choice y and v are not same. 
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For other interesting connections between graph and permutation limit theories refer to recent 
papers of Krai’ and Pikhurko [28], Glebov et al. [22, 23]. 


2.4. Empirical Degree Distribution of Random Permutation Graphs. Given a sequence 
of permutation graphs (G T „) n >i, the empirical degree distribution is 


1 n 

K (Cur„) ~ 


dn(i) ' 


i= 1 


Diaconis et al. [16] pointed out that if G nn —> W u , then n(G nn ) converges weakly to the distribution 
of the random variable 


WlQQ : = 


'[Oh] 2 


K{X, z)dv(z), 


where 2L = (Ad, Yi) is a random element in [0, l ] 2 with distribution z/, and K((xi,yi), (x 2 ,y 2 )) = 


l{(xi - x 2 ){yi - y 2 ) < 0 }. 

The degree distribution of a permutation graph is an immediate consequence of this discussion 
as explained below. 


Proposition 2.1. Let (tt,^^ be a sequence of permutations converging to a permuton v in 
([0,1] 2 ,^([0, l] 2 )), and G- Kn be the associated permutation graph. Then the empirical degree distri¬ 
bution 

1 n 

%X 1 +Y 1 -2F„(X 1 ,Y 1 ), (2.5) 

n f ur 

i=i 

where {X \, Y \) ~ v and F u is the distribution function of v. 

Proof. Let K : [0, l] 2 x [0, l] 2 [0,1] be A"((xi,yi), (x 2 ,y 2 )) = l{(aq - x 2 )(yi - y 2 ) < 0}. From 

the discussion before the proposition, n(G nn ) converges weakly to the law of 

Wi(*!,*!) = / 1 {(X 1 -x 2 )(Y 1 -y 2 ) < 0}du(x 2 ,y 2 ) = u([0,Ad] X \Y U 1]) + u ([0, H] x [X u l]). 

J [o,i ] 2 

The limiting degree distribution in (2.5) now follows from the fact that u has uniform marginals. □ 


3. Summary of the Results 

Let 7 T n E S n be any permutation. Define the random variable a n (i) := Ej=i ^ n (*,i), for 
i 6 [n]. Note that a n (i) E [0, i — 1] represents the number of edges in G Wn connecting the vertex 
i to vertices j E [i — 1]. The quantity a n (i) will be referred to as the back-degree of the vertex 
i. Similarly, one can define the front-degree of the vertex i as b n [i) := Jf l7 =i+1 D Wn (i,j). Note 
that b n (i) E [0, n — i\ and d n (i) = a n (i) + b n (i), where d n {i) is the degree of the vertex i. We are 
interested in the convergence of the degree process 

d n (t ) = d n {\nt])/n = a n (t ) + b n (t), (3.1) 

where a n (t) = a n (\nt\)/n and b n (t) = b n (\nt])/n, are the back-degree process and the front-degree 
process , respectively. 

Since [0, 1 ]( 0>1 ] is compact by Tychonoff’s theorem, a sequence of stochastic processes {Z n (t)} t£ ( 0 ,i] 
supported on [0,1] converges in law to a process {Z(f)} tg ( 0 q] (denoted by Z n {-) => Z(-)), if and only 
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if the finite dimensional distributions of Z n (-) converge to Z(-). A sequence of random measures 
p. n £ M. converges weakly in distribution to a random measure p, E M. (denoted by /j, n n), if 

Tn (/) := [ fdp n 4 [ fdn := //(/), (3.2) 

J[ 0,1] 2 J[ 0,1] 2 

for all continuous functions / : [0, l] 2 -A M. Finally, we denote by Jf(X), the distribution of the 
random variable X. 

Unlike the empirical degree distribution, the convergence of a permuton sequence (zAr n )n>i (de¬ 
fined in (2.3)), might not imply convergence of the degree process. An easy way to see this is to 
consider 7r n a sequence of deterministic permutations converging to A([0, l] 2 ), the Lebesgue mea¬ 
sure on the unit square, and a n a uniformly random permutation on S n . Define a new sequence of 
permutations K n which is n n if n is odd and a n if n is even. Clearly, v Kn converges to A([0, l] 2 ). 
However, for t E [0,1], the degree process d n {t) converges in distribution to the random variable 
Unif [t, 1 — t] along n even, and is a deterministic sequence of real numbers for n odd. 

The above discussion necessitates extra assumptions on the random permutations to ensure the 
convergence of the degree process. To this end, for -K n E S n and t E (0,1], define 

7 r n (f) := TT n (\nt])/n. (3.3) 

By this rescaling, 7r n : (0,1] i-a [0,1] is a stochastic process, hereafter referred to as the permutation 
process. The following theorem shows that the convergence of the permutation process implies that 
both the corresponding permutons and the degree process converge. To best of our knowledge, this 
connection between the convergence of the permutation process and the permutons is new, and 
might be of independent interest. In fact, under regularity conditions (discussed in Section 6) the 
two notions of convergence are equivalent. 

Theorem 3.1. Let n n E S n be a sequence of random permutations such that 

*»(■) ^ Z{-). (3.4) 

Then there exists a (random) measure p E Ai such that the permuton u 7Tn -a p, and the degree 
process 

d n (-)^D(-), (3.5) 

where 

D(t) =t + Z(t) — 2 Fpft, Z(t)), 
and F ^ is the distribution function of the measure p. 

The above theorem, which is proved in Section 4.1.1, will be used to determine the limiting degree 
process for various random permutations. Note that the limiting measure /i might be random. In 
this case, the finite dimensional distributions degree process can be dependent (see Example 4.1). 
However, for most of the examples considered in this paper the limiting measure in non-random, 
and the corresponding degree process has independent finite dimensional distributions. This is 
summarized in the following corollary and proved in Section 4.1.2. 

Corollary 3.2. Suppose the permutation process ir n (-) => Z(-), and the finite dimensional marginals 
of Z(-) are independent. 
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• Then converges to a non random measure fi G M. with law of ( X, Y) ~ fi as follows: 

X ~ Unif[0,1], and 2z?( Y\X = x) ~ jgf(Z(x)). 

• The finite dimensional marginals of the limiting degree process D(-) are independent. 


The above results imply that to determine the convergence of the degree process, it suffices to 
verify the convergence of the permutation process. This requires the permutation process to have 
some regularity, which is formalized in Section 6. The regularity conditions ensure that the finite 
dimensional distributions of the permutation process are “equicontinuous”, which can be verified 
easily for exponential models on permutation, like the Mallow’s model [32] and Spearman’s rank 
correlation models [13]. 


3.1. Uniform Random Permutation. A uniformly random permutation graph G nri is the per¬ 
mutation graph associated with a uniformly random permutation n n G S n . The degree distribution 
for a uniformly random permutation graph can derived easily from Theorem 3.1. To this end, de¬ 
note by [a, b\, for a,b G M, the interval [aA6,aV6], where aAb := min{a, b\ and aVh := max{a, b\. 
Moreover, Unif(a, b), for a, b G [0,1], denotes the uniform distribution over [a, b\. 


Corollary 3.3. Let ir n G S n be a uniform random permutation and d n (l), d n (2),... d n (ri) be the 
degree sequence of the associated permutation graph G nn . Then for indices 0 < r\ < r 2 < • • • < 

r s < 1, 


d n {\nr{\) d n (\nr 2 ~\) d, 




n n n 

where D\, D 2 , ■ ■ ■, D s are independent and Di ~ Unif(rj, 1 — r*), for every i 6 [s]. 


Proof. For a uniformly random permutation, the permutation process 7r n (-) Z(-), where Z(t) is 

independent Unif[0,1], for all t G [0,1]. Moreover, the permuton v Vri converges to the Lebsegue 
measure on [0, l] 2 . Therefore, for a G [0,1], by Theorem 3.1 

D a = a + U — 2 aU = (1 — a)U + a(l — U) ~ Unif(a, 1 — a), 
where U ~ Unif(0,1). □ 


Figure 1 shows the degree proportion of the labelled vertices in the permutation graph associated 
with a random permutation of length n = 10 ,r> . The symmetry in the figure around the mid-vertex , 
that is the vertex labeled n/2, is because the distribution of d n (i) and d n (n + 1 — i) are the same 
for every i G [n]. This is confirmed by the above Corollary, which shows that for r G [0,1], the 
degree proportion of the vertex labelled \nr~\ converges to the uniform distribution over the interval 
[r, 1 — r]. The shrinking length of this interval as r approaches 1/2, explains the fan-like structure of 
the degree proportion in Figure 1. Moreover, when r = 1/2 this implies that d n (\n/2])/n converges 
in probability to 1/2, as can be seen from Figure 1. 

It follows from Corollary 3.3 that the degree proportion of the mid-vertex, that is, the vertex 
labelled [n/2], converges to 1/2 in probability. Therefore, it is reasonable to expect a central limit 
theorem for d n ([n/2])/n around 1/2 after an appropriate rescaling. This is detailed in the following 
theorem and illustrated in Figure 2(a). The proof is given in Section 5.3. 
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Sampling Distribution 


Empirical Degree Distribution 



(a) 


(b) 


FIGURE 2. (a) The sampling distribution of the degree of the mid-vertex (the vertex labeled 
n/2), for the permutation graph associated with a random permutation of length n = 10 4 , 
repeated over 10 5 samples, (b) The empirical degree distribution of the permutation graph 
associated with a random permutation of length n = 10 5 . The blue curve represents the 
density of the limiting distribution f(z) = — log \ 2z — 1|. 


Theorem 3.4. Let Ti n £ S n be a uniform random permutation and d n (l), d n (2),... d n {n) be the 
degree sequence of the associated permutation graph G nn . Then 

yfr (MY 2 !! - 0 A N( 0,17(1 - 17)), (3.6) 

where U ~ Unif[0,1]. 

The empirical degree distribution of a uniformly random permutation graph is a direct conse¬ 
quence of the Proposition 2.1. The limiting density of the empirical degree distribution is depicted 
in Figure 2(b). The density is supported on [0,1] and has an interesting shape: it vanishes at the 
end points, and blows up to infinity at £ = 1/2. 

Corollary 3.5. Let ir n £ S n be a uniformly random permutation and G- Kn the associated permuta¬ 
tion graph. Then the empirical degree distribution 

1 n 

«(G 7r J := - V 6 dn(i) 4 Z := (1 - U)V + U(l - V ), (3.7) 

n Z —' ^nT~ 

1=1 

where U. V are independent Unif(0, 1). Equivalently, Z has the same distribution as Unif(V, 1 — V), 
where V ~ Unif(0,1), and has a density with respect to Lebesgue measure given by fz(z) = — log |1 — 
2z\, for 0 < z < 1. 

The degree distribution of a uniformly random permutation graph has different variability de¬ 
pending on the label of the vertex. In fact, the next corollary shows that a uniformly random 
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permutation graph is rather irregular: a typical vertex has degree around nj 2 , however the mini¬ 
mum and the maximum degrees are far apart: 

Corollary 3.6. For a uniformly random permutation ir n , n~ i S(G 7Tri ) —> 0 and n -1 A ((?„■„) —> 1, 
where A(G nn ) and S(G irn ) are the maximum and the minimum degree of G- Kn , respectively. 

Proof. To prove the result for the minimum degree it suffices to show that for any 0 < e < 1/2, 
lim^oc P(5(G 7rn ) > ne) = 0. To this end, choose r > 1 arbitrary, and real numbers 0 < t\ < t 2 < 
■■■ <t r < e/4. By Corollary 3.3 

lim P n f min d n (\nti}) > ne) = TT V ° /' < (^ 
iho o \l<i<r ) 1-2 ti \l-e/2J 

Since this is true for every r and for all e > 0 fixed, letting r —> oo gives the result. The result for 
the maximum degree can be proved similarly. □ 

In light of the above corollary, it is natural to ask whether 5(G nn ) converges to a non degenerate 
distribution after appropriate rescaling. The following theorem, which is proved in Section 8 shows 
that the minimum degree of a uniformly random permutation graph is a Rayleigh distribution. 


Theorem 3.7. For a uniformly random permutation ir n , 

*(GQ 3 r 

sfn 

where T is the Rayleigh distribution with parameter ^, that is, P(T > 7 ) = e -0 " 2 for all 7 > 0. 

3.2. Mallows Random Permutation. One of the most popular non-uniform model on permu¬ 
tations, which has applications in statistics [32], is the Mallows measure. For /3 G M, denote by 
7 r n ~ the Mallows random permutation over S n with probability mass function 

_o A (°~) 
e r* n 

m P,n(. a ) := zp 

T,aes n e n 

where A(<r) = |{(z, j) : (i — j)(Tr n (i) — vr n(j)) < 0}| is the number of inversions of the permutation a. 
The uniform random permutation corresponds to the case /3 = 0. It is easy to check that for 717 e S n 
chosen uniformly at random, v nn converges weakly in probability to Unif(0,1). This was generalized 
to Mallow random permutations by Starr [40]. Using this and Theorem 3.1 we compute the limiting 
density of the degree proportion in a Mallows random permutation. The limiting density exhibits 
interesting phase transitions depending on the value of f3. This is summarized in the following 
theorem and proved later in Section 7. 


Theorem 3.8. Let (3 G M and ir n ~ M^ n . If d n (l),d n (2),... d n [ri) is the degree sequence of the 
associated permutation graph G 7Tn , then the degree process 

dn(-) ^ Z(-). 


More precisely, for indices 0 < 77 < r 2 < 
'd n (\nr 1 ]) d n ([~nr 2 ]) 


<r s < 1 , 
d n {\nr s ]) 


n 


n 


n 


^ ^ (Di j( a, T>2,/3, ■ ■ ■, D Sj p), 
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FIGURE 3. Density of D a ^ for /? = 2 and (a) a = 0.1, and (b) a = 0.55. For /3 = 2, a c ((3) = 
0.28311. Since 0.1 ^ (0.28311,0.71689), the density (a) of -D 0 .i ,2 is a continuous function 
supported on [0.1,0.9]. On the other hand, 0.55 € (0.28311,0.71689), and the density (b) 
of Do. 55,2 is a piecewise continuous function with the pieces supported on [0.38,0.45) and 
[0.45, 0.55], with a discontinuity at 1 — a = 0.45. 


where D i ;i g, D- 2 ,p ,..., D s $ are independent, and 
1: if a ^ [a c (/3), 1 - a c (/3)], /ias density: 

p e \p{ a ~ z ) 

9a,p (z) = -r, ze[a,l-a\] 

(1 - e~P) (e^ a+z ) - ePR(a, /3)) 2 

2: if a 6 (a c (/3), 1 — a c (/3)), D aj/ g /ias density: 


9a,p(z) 


where 


' p e \p( a ~ z ) 

(1 -e-^)(e' 3 (“+ z )-e/ 3 it(o,/3))2 
2/3 e hP( a ~ z '> 

(l-e~P)(eP ( . a + z '>-ePR(a,P))^ 


for z E [a, 1 — a], 


D(a,/3) 


4 (e^ 3 - e a/3 ) (e a ^ - l) 
(e / 3 - l ) 2 


and a c ((3) 


1 logcosh(/3/2) 

2 


(3.8) 


(3.9) 


(3.10) 


The above theorem gives the limiting distribution of the degree process of the permutation graph 
G w „ associated with a Mallows random permutation 7r n ~ Mp n . For ft 6 M fixed, the limiting 
distribution of d n (\na])/n has a phase transition depending on the value of a E [0,1]. There exists 
two critical points a c (f3) and 1 — a c (/3), such that for a [a c (/3), 1 — a c (/3], the limiting density 
of D a ,p is a continuous function supported on [a, 1 — a]. However, for a in the critical interval 
(a c (/3), 1 — a c (/3)), the density of D a>l g breaks into two piecewise continuous parts on the intervals 


1 — a + — log R(a, f3 ), a A 1 — a 


and (a A 1 — a, a V 1 — a], 


















DEGREE SEQUENCE OF RANDOM PERMUTATION GRAPHS 


13 



FIGURE 4. The support of D aA 3 , for /? = 2 and a £ [0,1]: The points of tangency of the 
curve b = 1 — a + 4 log R(a, (3) with the two straight lines b = a and b = 1 — a are colored 
green. The corresponding red points on the i-axis are the critical points a c (/3 ) and 1 —a c (/3). 
For a fixed a £ [0,1], the support of D a j 3 is the interval intercepted by the vertical line at a 
either between the two straight lines (if a ( a c {/3 ), 1 — a c (/3 )), for example, when a = ai), 
or between the curve and one of the two straight lines (if a £ (a c (/3), 1 — a c (f 3)), for example, 
when a = 02 ). 


with a discontinuity at the point a A 1 — a. Plots of the limiting density of D a ^ are shown in Figure 
3 for /3 = 2 and a = 0.1 and a = 0.55. The changes in the support of D a ^ for values of a in the 
critical interval is depicted in Figure 4. 

The critical curves f3 eA a c (/3) and (3 eA 1 — a c (/3) are shown in Figure 5. For a fixed £ M the 
critical interval (a c (/?o), 1 — a c (/3o )) is the interval between the 2 curves intercepted by the vertical 
line at (3$. Note that for (3 = 0, the a c ((3 ) = 1 — a c ({3) = 1/2, that is, the critical interval is 
empty. Therefore, for a uniform random permutation, the limiting density has no phase transition, 
as elaborated in Corollary 3.3. 

The phase transition in the density of D a p can be reinterpreted by fixing a £ [0,1] and varying (3\ 
Theorem 3.8 shows that for a £ [0,1] fixed, there exists a critical point (3 c (a ) (obtained by solving for 
/3 in a c ((3) = a) such that for (3 £ [0,/3 c (a)], the density of D a ^ is a continuous function supported 
on [a, 1 — a]. However, for f3 > (3 c {a ) the density of D a ^ breaks into two piecewise continuous parts 
with a discontinuity at the point a A 1 — a. If (3 = 1/T denotes the inverse temperature, then this 
phenomenon is the effect of replica symmetry breaking in statistical physics as one moves from the 
high temperature to the low temperature regime. 

4. Degree Distribution of Random Permutations 

In this section we prove the convergence of the degree process whenever the permutation process 
converges. Recall that the Kolmogorov-Smirnov distance on the space of probability measures on 
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FIGURE 5. The critical transition curves a c (f 3) and 1 — a c (/3 ) for f3 £ [—10,10]. For a 
fixed P 0 £ K the critical interval (a c (/l 0 ) ) l — a c (/3 0 )) is the interval between the 2 curves 
intercepted by the vertical line at /3q. 


[ 0 , l ] 2 is defined by 

\W - v\\ks ■= sup \F u (x,y) - F^{x,y)\, 

0<x,y<l 

where F v and Ff are the bivariate distribution functions of v and g respectively. Note that con¬ 
vergence in Kolmogorov-Smirnov distance implies weak convergence, but not conversely. Finally, 
for ir n £ S n define the empirical permutation measure 

^ 7r n ~ — ^ 1 ™n(i) 

ie[n] Kn ’ n 

It is easy to check that the permuton u nn associated with the permutation 7 r n satisfies 11 u 7rn — 
v-K n 11ks 0, and so a sequence of permutations 7r n £ S n converges to a permuton v if and only if 
z>jr n converges weakly in probability to v. 

The following theorem shows that the convergence of the permutation process implies the con¬ 
vergence of the permutons. This connection between the two types of convergence might give new 
insights into the permutation limit theory described in Section 2. 



Theorem 4.1. Let n n £ S n be a sequence of random permutations such that 

*»(■) ^ Z(-). (4.2) 

Then (for n5 foi( - )) converges jointly weakly in distribution. In particular, there exists a (possibly 
random) measure \x £ Ai. such that the permuton sequence (for„)ra>i converges in distribution to 
H £ At. 

Proof. Fix a continuous function / : [0, l] 2 -» [0, 1], positive integers a,b > 1, continuous functions 
gi,g 2 , - ■ ■ ,gb ■ [0,l]i->K, and real numbers si, S 2 , ■ ■ ■ ,Sb£ (0, 1]. The joint a-th moment of v n . n (/) 
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and {gi(TT n (si)),g 2 (Tr n (s 2 )), ■ ■ ■ ,g n (^n(sb))) is 

5 ' 1 


E [KADUdjMsj))] = —E ^ n f 


3 =1 


3 =1 


ij (J'j ) 


n n 


Y[9aMsj)) 


3 = 1 


( b « a 

Ylsj^nisj)) / TT f{x j ,Tr n {x j ))dx j ] +o(l), (4.3) 
7=i 


y? =I 

where the last equality follows from the uniform continuity of /. 
Now, by assumption (4.2), 

b a b a 

n 9j{*n(sj)) n f{xj, ir n (xj)) 4- Yl 
3=1 3=1 3=1 3=1 

Therefore, by the Dominated Convergence Theorem 
b 


j™o e n®M-i))n/( 


•Ej 5 '^nyKj 




U = 1 


3=1 


D=! 


3=1 


The RHS above is measurable in (aq, X 2 , ■ ■ ■ , x a ), as it is the limit of measurable functions. Another 
application of Dominated Convergence Theorem gives 

e i n«w-i))n/(-i.^)) l n dxj + o(l), 


E ( (/) II 9jMSj)) 

3=1 


'[ 0 , 1 ]“ 


U = 1 


3=1 

a 


3=1 


•'[°, 1 ]“ L =1 j=l / 3=1 


Since (4.4) holds for all choices of / and gi,--- ,gb, the joint convergence in law of (z/^,vr n (-)) 
follows. Since the space Ad is closed, implies P^ n —> g, for some random measure fj, G M.. 

Finally, since \\v-m ~ fv n ||i<:s —> 0 in probability, (&v„)n>i converges in distribution to the per- 
muton n G M. □ 

4.1. Proofs of Theorem 3.1 and Corollary 3.2. The above theorem can now be used to prove 
the convergence of the degree process. 


4.1.1. Proof of Theorem 3.1. By Theorem 4.1 we have 

(£*„,*■ n(0) 

where ju G M is a random measure such that u nn —> n, and Z{.) is a stochastic process on (0,1] 
such that 7r n (-) —> Z(-). Thus, fixing k > 1 and 0 < t\ < t2 < ■ ■ ■ < t k < 1 we have 

,7Tn(4)) (n,Z(tl),--- ,Z(t k )). 

Applying Skorohod’s representation theorem on the separable metric space M x [0, l] fc (see Billings¬ 
ley [4, Theorem 6.7]), without loss of generality assume that the above convergence happens almost 
surely. 
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Now for any t G (0,1], 

1 \nt] 

a n (t ) = -^l{7 T n (o) > TT n (\nt])} = iv„([0, I" nt] /n] x (7r n (f),l]) 

a =1 


= ^([0,t] X (7T n (f),l]^) +0(1) 
=t ~ F Vn n (t,TT n {t)) + o(l), 

(4.5) 


where the last step uses \\v^ n — iAi-JIks = o(l). By a similar argument, 

&n(0 = t r„(t) - (f, 7T n (t)) + o(l). (4.6) 

Combining (4.5) and (4.6), for any 1 < i < fc we have 


dnifi ) — O^nifi) 4“ b n (ti) —4“ TTn(t'i) 2 F VlT ^{ti^TT n {tff) 4“ o(l) 

— ti 4“ 7T n (tj) 2F^(ti, 7T n (tj)) 4“ 11 ^TT n mIIks 4~ o(l) 

4t i + Z(t i )-2 J F / ,(i i ,Z(t i )), 


where the last step uses ||z/^ n — /z||ks 4 0 (see Hoppen et al. [26, Lemma 2.1]), and the fact that 
the function is continuous in each co-ordinate when the other co ordinate is held fixed. Indeed, 
this follows from the observation that any /i£ A4 has continuous marginals. Thus we have 

{d n (h),-- ■ ,d n {t k )) 4 (ti 4- Z(ti) - 2Fn(t\, Z(ti)), ■ ■ ■ ,t k + Z(t k ) - 2F fJ (t k , Z(t k ))), 
from which finite dimensional convergence of d n (-) follows. 


4.1.2. Proof of Corollary 3.2. From (4.3) and (4.4), it follows that for any continuous function 
/ : [0, l] 2 -A [0,1] and a > 1, 

lim EzA r „(/) a = f E 

n^oo 7 [ 0 , 1 ]“ 

Under the assumption of independence, the RHS above equals ^ E f(x, Z(x))dx j , which implies 

that u Wn (/) converges in probability to the non random quantity 

Ef(x,Z(x))dx = E fi f(X,Y), 

where (X , Y) ~ ^ is as required. Since this holds for all continuous functions / the desired 
conclusion follows. 

The independence of the finite dimensional marginals of Z(-), implies the same for the degree 
process D(-) by (3.5). 


a \ a 

n f( x i,z( x j)) n dxj- 

J =i / i =1 


4.2. A Dependent Degree Process. Even though theorem 3.1 allows for n to be random, in 
most examples in this paper /i turns out to be non-random and the corresponding degree process 
has independent finite dimensional distributions. In this section we construct a sequence of random 
permutations where the limiting permuton is random and the finite dimensional distributions of 
the degree process are not independent: 
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Example 4.1. Suppose W n is a uniform random variable on [n] := {1,2, and 7 r n G S n 

defined by 

n n (i) := (i + W n - 1 mod n) + 1. (4.7) 

Note that 7 r n is a cyclic shift of the identity permutation, where the length of the shift is chosen 
uniformly random. 

Proposition 4.2. Let ( 7 r n ) n >i be a sequence of random permutations as defined in (4-V- Then 
the degree process 

d n (t) ^ D(t) := W ■ 1 {W + t < 1} + (1 — W) ■ 1{W + t > 1}, (4.8) 


where W ~ Unif[0,1]. 

Proof. We will first show that the permutation process 7 r n (-) converges weakly in distribution. For 
s > 1 and let gi,g 2 , m ■ ■ ,9s be continuous functions on [0,1]. Then, 

/ ^ x ri 


%i^n(ii)J • • • 9s[TTn{is)) 


gi(ii + u mod 1 ) • • • g s {i s + u modi )du. 


Hence, 7 r n (-) => Z(-), where Z(-) is a stochastic process defined by 


Z(t ) = W + t mod 1, with W ~ Unif[0,1]. 


(4.9) 


By Theorem 3.1 and by (4.3) we have u nn —> fa such that for any continuous function / on the 
unit square, 

„ a a 

d x j = If f( x ji w + x j mod 1) T | dxjdu 

J ^ a+1 jA 

1 / /•> \ a 
/ f(x,x + u mod l)dx ) du: 

Therefore, the limiting measure fi is random and has the following law: For every s G [0,1], let n s 
be the joint law of (V, s + V mod 1), where V ~ Unif[0,1]. Then g = kw with W ~ [0,1]. 

To compute the limit of the degree process, we compute the distribution function of n s . In this 
case, with U ~ Unif[0,1] and 0 < a, b < 1 



mfr = 




E | l\f( Xj ,Z(xj)) 
1 



F Ks (a, b ) =P(f7 < a, U + s mod 1 < b) 

=P (U < a,U + s < b) + F(U < a, l<U + s<b + l) 

= min(a, b — s)+ + min(a + s — 1, 6)+, (4-10) 


which implies 

F Kw (t,Z(t))=mm(t,(W + t mod 1) — W)+ + min(t + W — 1, W + t mod 1)+ 

=t ■ 1 {W + t<l} + (t + W- 1)1 {W + t > 1}. (4.11) 

Therefore, by Theorem (3.1) d n {-) D(-) where: W ~ U[ 0,1], and 

D(t) :=t + (W + t mod 1) — 2 F Kw (t, W + t mod 1) 

=W-l{W + t< 1} + (1 — W) -\{W + t> 1}. 


(4.12) 
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In this case the finite dimensional distributions of the permutation process (4.9) and the degree 
process (4.12) are not independent, and the limiting permuton is random. □ 


5. Uniformly Random Permutation Graph 

Limiting properties of permutation statistics associated with a uniformly random permutation 
are widely studied. The number of edges \E(G nn ) \ is the number of inversions in the permutation 7 r n . 
For a uniformly random permutation the distribution of the number of inversions is Xa=i where 
the random variables X % are independent and uniformly distributed over the set {0,1 ,i — 1}, 
for every i 6 {1,2,... ,n}. Normal approximations to \E(G 77n )\ can be proved using these results 
and standard versions of the central limit theorem (refer to Fulman [21] for a proof using the 
method of exchangeable pairs). The largest clique in a permutation graph corresponds to the 
longest decreasing subsequence in the permutation. Similarly, an increasing subsequence in a 
permutation corresponds to an independent set of the same size in the corresponding permutation 
graph. Asymptotics for the maximum clique and the independent set in G nn follow from the seminal 
work of Baik et al. [3] on the length of the longest increasing subsequence in a uniformly random 
permutation. 

The convergence of the degree process of a uniformly random permutation graph (Corollary 3.3) 
is a direct consequence of the general theorem. The empirical degree distribution can also be easily 
derived. 


5.1. Proof of Corollary 3.5. When 7 r n is a uniformly random permutation, the limiting permuton 
is v = Unif(0,1) x Unif(0,1) and the result follows from Proposition 2.1 by direct substitution. 

To get the density of the limiting random variable, let Z := (1 — U)V + U{1 — V) where U, V 
are independent Unif[0,1]. For z < 1/2, conditioning on U the distribution function of Z can be 
calculated as 




z — (l — u) 
2 u- 1 


du. 


1 - 2 u 

Simplifying and differentiating the above expression with respect to z gives the desired density for 
z < 1/2. For z > 1/2, the density can be derived similarly. The density vanishes at the end points, 
and blows up to infinity at z = 1 / 2 . 


5.2. A Hypergeometric Estimate. For studying the degree sequence of a uniformly random 
permutation graph, properties of the random variables a n (i ) and b n (i) (defined in Section 3) will be 
needed. To this end, recall the hypergeometric distribution: A non negative integer valued random 
variable X is said to follow the hypergeometric distribution with parameters ( N , M, r) if 


P(X = x) 



for x E [max{0, r + M — N}, min{M, r}] 


where N > max{M, r}. 

Note that (i— 1) — a n (i)+b n (i) = n n (i) — 1, thus giving the simple relation b n (i)—a n (i) = n n (i) — i. 
Using this relation, the following proposition gives a concentration result for d n (i) around its 
conditional mean given TT n (i). 
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Proposition 5.1. The conditional distribution of a n (i)\{n n (i) = j} is hypergeometric with param¬ 
eters (n — 1 ,i — 1 ,n — j). Consequently, for R > 0 


d n (i) - 


(i - 1 )(n - 7r„(z)) - (7r n (i 

- l)(n — z) 

n — 1 



> R 


TTnW = j < 2 e 2 n 


Proof. Given ir n (i) = j and a n (i) = a, b n (i) = a + j — i =: b, and so 

(i - l)\(j - l)\(n - i)\(n - j)\ 0 a 1 ) ( ( (n-j)-a } ) 


P(a n (z) = a\n n (i) = j ) = 


a!(z — 1 — a)! 6 !(n — i — b)\(n — 1 )! 


( n_1 ) 

\n—iJ 


and so a n (i) follows the hypergeometric distribution with aforementioned parameters. 

Therefore, E(d n (z)| 7 r n (z) = j ) = E(a n (z)) + E (b n (i)) = THvl+L 1±L'~ 1 >( rt ~l) . To prove the second 

conclusion note that 

(* - 1 )(n-j) + (j - l)(n — i) 


d n {i) 


n — 1 


> R 


a n (i) 


(*- 1 ){n- j) 


n — 1 


i? 

> 2 - 


An application of the bound in [39] now gives the desired conclusion. 

5.3. CLT for the Mid-Vertex: Proof of Theorem 3.4. Let Z n = y/n 

fixing 5 > 0 

1 n 

nz n < X) = - V'p(z n < x\n n (\n/2~\) = j) 

n zJ 

3 = 1 

= - E P(-^n < x\^ni\n/2\) = j) + 0(5) 

71 t 


d n (\n/2\) _ 1 
n 2 


□ 

Now, 


tt,< 5 <n( 1—<5) 


- E 


'a n (\n/2\) - 


(K21-1 )(n-j) 

n— 1 


n<5<j><n(l—(5) 


l 

77- 2 


< An(x, j) 7 T n ([n/ 2 ]) = jj +0(5), 

(5.1) 


where A n (x,j) satisfies linin^oo ma ’x n S<j< n (i-S) |A n (aj,j) — x/2| = 0, for all 
By Proposition 5.1, 

rt(\n/2-}),j) := Var(o Tl (fn/ 2 ])| 7 r n (["n/ 2 ])) = j) 

(\n/2\)~ 1 )(j - 1 )(n — \n/2]))(n-j) ^ 

= - 7 -TVs?-R 7 -> C(5)n, (5.2) 

(n — l) z (n — 2 ) 

for some C(5) > 0, and for all j E [n5,n( 1 — <5)]. Using the Berry-Esseen theorem for hy¬ 
pergeometric distribution [29, Theorem 2.2], there exists a universal constant C such that with 
C’(5) := C/y/W) < oo, 


a n (\n/2\) - 


(r»/ 2 l — l)(n— j) 
n— 1 


l 

n 2 


< A n (x, j) TT n (i ) = j - $ V 7 ™ 


An(^,j) 

o- n (fn/ 2 l, j) 


<m (M) 


— l 

n 2 


Finally, note that 


max 

nS<j<n(l—<8) 


°l(W/2\,j) j(n — j) 


n 


4 n 2 


= o(l), 
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where the o(l) term goes to zero as n -> oo. Moreover, since the function is uniformly continuous 
on M, 


th I /— j ) 

max $ y/n ■ —————— 

nS<j<n(l—S) \ a n (\n/2\,j) 

Combining (5.1), (5.3) and (5.4) we have 

1 


- 


VtiMi 1 ~ i/ n ) 


= 0 ( 1 ). 


(5.4) 


P (Z n <x) = - V $ 
n 


\/{jln){l- j/n) 


nS<j<n(l—5) 

On taking limits as n —> oo followed by 5 —> 0 we have 

lim P (Z n < x) = f ( — X 

n ^°° Jo \ -y/u(l — u) 

which completes the proof of the theorem. 


+ o(l) + 0(6). 


du, 


6 . Convergence of the Permutation Process 

The convergence of the permutation process requires some regularity assumptions. In this section 
we introduce the notion of equicontinuity for a sequence of random permutations, and verify this 
for most standard exponential models on permutations. 

6.1. Equicontinuous Permutations. We begin by recalling few definitions: 

Definition 6.1. Let & be a family of functions from [0,1] to M. The family & is equicontinuous 
at a point xq G [ 0 , 1 ] if for every e > 0 , there exists a <5 > 0 such that 

| f (.To) — /(a;)| < £ for all / G & and all x such that \x — xo| < 6. 

The family is equicontinuous if it is equicontinuous at each point of [0,1]. 

Denote by C[0,1] the set of all continuous functions from [0,1] to M. Similar to the notion 
of equicontinuity of a class of functions, we introduce the following notion of equicontinuity of 
permutations. 

Definition 6.2. Let a > 1 and g\, g 2 , ■ ■ ■, g a £ C[0,1]. For any sequence 7 r n G S n of random 
permutations, define the function G gi : [0, l] a —> M 

G SSs 2 ,...,ff«( S l’ S 2 ,---,Sa) =E ^^(^(Si))^ ■ ( 6 . 1 ) 

A sequence ( 7 r n ) n >i of random permutations is said to be equicontinuous if the family {Ggl] g2> ... t g a } n >i 
is equicontinuous for all 91 , 52 ,... , g 0 G C[0, 1] and a > 1. 

Proposition 6.1. Let Tr n G S n be a sequence of random equicontinuous permutations such that 
v TTn L- Then the permutation process 


*n(0 ^ Z{.) 
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where the finite dimensional distribution of Z(-) is as follows: Let (Xi, Yi), (X 2 , Y 2 ), ■ ■ ■ ( X a , Y a ) be 
independent draws from the random measure n E Xi. Then 

J?(Z( Sl ), Z(s 2 ), • • • , Z(s a )) ~ JS?(Yi|Xi = ai, y 2 |X 2 = a 2 , ■ ■ ■ , Ya\X a = So), ( 6 . 2 ) 

for 0 < si < s 2 ■ ■ ■ < s a < 1. 


Proof. Fix a > 1. For every n > 1 define a collection of random variables {{Uj <n , Vj,n)}j=i> where 
{Uj,n}j = 1 are i.i.d. Unif[0,1] independent of 7 r n , and V), n := iT n (Uj^ n ). Let / 1 , • • • , fa ■ [0, l] 2 —> M 
be continuous functions. Now, as in (4.3) 



a \ a 

Y\fj{Xj,TT n ( x j)) J \dXj 

4=1 J 


7=1 


(6.3) 


Since z/^ —> /a, the RHS of (6.3) converges to 

CL n CL 

/ T[fj(x j ,y j )dfj,{x j ,y j ). 

t=i 

Thus, the joint law of {(Uj jTl , Vj, n )}“ =1 converges to {{Uj, Vj)}“ =1 , where {Uj}j =l are i.i.d. Unif[0,1] 
independent of /i, and given both {Uj = Uj}j= 1 and g,we have { Vj}“ =] are mutually independent 
with V) M(>j). 

Observe that for g\. g 2 , .... g a £ C[0,1], 


E 



£/l,n 


— ®i > ■ j LJ a n 



=En«w-i)) 

7=1 


(6.4) 


and 

) j , (6-5) 

To show that 7r n (-) =4> Z(-), it suffices to show that the RHS of (6.4) converges to the of RHS (6.5) 
are equal. Therefore, it suffices to prove that the corresponding conditional distributions converge: 


E 


n»<u> 

\i = 1 


U\ — ^ 1 ? * 5 U a — s a I — IE 



3 = 1 


3 \ °J 


2mn}U\{Uhn}U ) ^ 

This follows by the equicontinuity of (7r n ) n >i and an application of [ , Theorem 4], since [0,1]“ is 

compact. □ 


The following proposition gives a criterion for equicontinuity of permutations. This can be easily 
checked for many non-uniform models on permutations. To this end, let x = (xx,x 2 ,--- i x s) £ 
(0, l] s and define 

^(*i,*2, • ■■is,x) := {ir n £ S n : n n (\nx{\) = h,-- • ,7T n (f ~nx s ]) = * s }. 


( 6 . 6 ) 
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Proposition 6.2. Let n n G S n be a sequence of random permutations. Suppose for all s > 1 and 
1 < i\ < ii < ■ ■ ■ < i s < n, the following holds: 


lim lim sup 

5—>-0 n—^oo (x,y)eB(5) 


P(7T n G n(ii,i 2l - ■ -i s ,x)) 
P(7T n G • ■■is,y)) 


(6.7) 


where B s {5 ) = {a^i, - - - , x s ,yi,-- - , t/ s G (0,1] : |x s — y s | < 5}. XTien f/ie permutation sequence 
{^n)n >l equicontinuous. 


Proof. Fix e > 0. Let x = (xi, X 2 , ■ ■ ■ , x s ) G (0, l] s and y = (yi, y 2 , ■ ■ ■ , y s ) G (0, l] s . Denote by 
K n _ s the collection of s-tuples ■ ■ ■ ,i s which are all distinct. 

To show equicontinuity, fix s > 1 and g±, g 2 , ■ ■ ■, g s G C[ 0, 1]. Let M := sup 3 . 1 ... a . sg ( 0 ,i] F[/=i 3j( x j )■ 
Now, using (6.7) 

SU P - G( £ ] , 92 ,..., 9 s(y)\ 

x,y£B 3 (6) 


< 


sup 

x,yeM5) {iui2 ... tis)eKny 


S 


rib 


< Me, 



• |P(7T n G n(ii,*2,--- ,is,x)) -P(vr n G fi(n, *2, - - - ,*s,Z/))| 


( 6 . 8 ) 


for large enough and appropriately chosen 5. This implies that the sequence of permutations 
{^n)n >l is equicontinuous. □ 


6.2. Exponential Models on Permutations. Let 9 G M, and T n : S n —> M be any function. 
Suppose 7 T n G 5„ is a sequence of random permutations with probability mass function 


e 9T n (cr) 


y e eT n {a)' 


(6.9) 


One of the most common exponential models on permutations, is the Mallows model, where T n (a ) 
is the number of inversions of a scaled by n. The limiting degree process of the Mallows random 
permutation is explicitly computed in Section 7. Here, we determine a simple criterion for the 
convergence of the degree process of a sequence of random permutations distributed as (6.9). To 
this end, consider the following technical definition: 


Definition 6.3. For any two fixed vectors x = (x\,X 2 , • • • , x s ) G (0, l] s and y = (t/i, t/ 2 > ■ ■ ■ , y s ) G 
(0, l] s , define the bijection 


®x,y ■ fi(*l,*2, ■ ■ -is,x) ->• ■ -is,y) 


as follows: for each ir n G D(ii, 12 , ■ ■ ■ i s ,x), define its image n n as 


TT n {k) 


ij = Tr n (\nxj~\) 
< n n (\nyj] 
TT n (k) 


if k = | \nyf\ for j G [s] 
if k = \nxj\ for j G [s] 
otherwise. 


( 6 . 10 ) 


Informally, <§>x,y takes a permutation n n and interchanges the [nxi], [nx' 2 ] ,... \nx s ~\ coordinates to 
the \nyi\, \ny 2 \, ■ ■ ■ \ny s ~\ coordinates to get 7f n . Using this bijection it is easy to get a sufficient 
condition for the convergence of the permutation process for exponential models. 
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Corollary 6.3. A sequence of random permutations (7r n ) n >i from (6.9) is equicontinuous whenever 

lim lim sup sup |T n (7r n ) - T n {^ x ,y(Fn))\ = 0, (6.11) 

<5-5.0 n-too ( x,y)eB(6 ) 7 r n eO(ii,... ,i s ,oc) 

where & x ,y is the bijection defined in (6.10). 

Proof. Fix e > 0. Let n n be the image of n n defined by the bijection Q x ,y hi (6.10). Then using 

( 6 . 11 ) 

P(7r n 6 Q(ii,i 2 , • • • ,i s ,x)) _ 1 

P(7T ri e n(ii,*2, • • • ,is,y)) 

for n large enough and appropriately chosen 6. Equicontinuity follows from Proposition 6.2. □ 


S7r n eQ(ii,i2,-.. ,i s ,x) QnfiiFn) 


Consider the following general class of 1-parameter exponential family on the space of permuta¬ 
tions S n , with probability mass function 







E 


cr&Sn 


«£ Uf 




( 6 . 12 ) 


where / : [0, l] 2 -A [0,1] is any continuous function. This is a special case of the model (6.9) 
with T n (a) = Y17= i f {n‘ • Popular choices of the function / includes the Spearman’s Rank 

Correlation Model : f(x,y) = —(x — y ) 2 and the Spearman’s Footrule Model f(x,y ) = —\x — y\. 
These models find applications in statistics for analyzing ranked data [13, 1 I]. Feigin and Cohen 
[20] gave a nice application of such models for analyzing agreement between several judges in a 
contest. For other choices of / and their various properties refer to Diaconis [ ]. Consistent 

estimation of parameters in such models has been studied recently by Mukherjee [34]. 

Using Corollary 6.3 it can easily shown that any sequence of random permutations (7r n ) n >i 
distributed as (6.12) is equicontinuous, that is, their corresponding degree process converges: 


Corollary 6.4. Fix 0 € M and f : [0, l] 2 [0,1] be any continuous function. Let ir e S n be a 

sequence of random permutations distributed as (6.12). Then (7r n ) n >i is equicontinuous and the 
degree process d n (-) converges. 


Proof. The convergence of the degree process follows from the equicontinuity (Proposition 6.1). To 
prove equicontinuity, let T n (a) = /(*/ n > °'(*)/ n )- Fi x 6 > 0 and let (x,y) G B(6). Using the 

bijection (6.10) we get 


S 

\T n {7r n ) — T n (fr n )\ <2^ 

3 = 1 

<2 s sup \f(xuy) - f(x 2 ,y)\, 

\xi-x 2 \<S+f, 

2/S [0,1] 

which goes to 0, after taking n —> oo and <5 —A 0, by the continuity of /. 


/ 


\nxj 1 n n (\nxj~\) 


n 


n 


~f 


\nyj] 7r n (\nxj]) 


n 


n 


(6.13) 

□ 
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7. Mallows Random Permutation 


For /3 G M and n > 1, the Mallows measure over permutations 7r n G <S n is given by the probability 
mass function 

mp tn ((j) : = 


= " n. 


-P 




where A(cr) = |{(z, j) : (i — j)i^n 1 (i) ~ 7r n 1 U )) < 0}| is the number of inversions of the permutation 
1 7. Diaconis and Ram [17] studied a Markov chain on S n for which the Mallows model gives the 
limiting distribution. Tail bounds for the displacement of an element in a Mallows permutation 
was studied by Braverman and Mossel [ ]. Recently, Mueller and Starr [ ] and later Bhatnagar 

and Peled [6] studied the length of the longest increasing subsequence in a Mallows permutation. 
Consistent estimation of parameters in exponential models on permutations has been studied by 
Mukherjee [34]. 

The uniform random permutation corresponds to the case fd = 0. It is well known that for 
ir n G S n chosen uniformly at random, v nn converges weakly in probability to Unif(0,1). This was 
generalized to Mallow random permutations by Starr [40] . 

Theorem 7.1 (Starr [ 10]). Let ir n G S n be a Mallows random permutation with parameter f3. Then 
the empirical permutation measure D 7T - n converges weakly in probability to a random variable which 
has density in [0, l] 2 given by 

(5 sinh f | 


mp(x,y) := - 

2 ( ex P (f) cosh (^( x _ y)) ~ ex P (-?) cosh + y ~!))) 

and distribution function 

2exp (—P/3(a + 6-1)) (sinh sinh 


2 ’ 


(7.1) 


Mp(a, b) = - ^ log 


sinh ( f 


- 1 


(7.2) 


Proof. The proof of (7.1) can be found in Starr [40]. The expression for the distribution function 
(7.2) follows by directly integrating the density mp in (7.1). □ 

By the above theorem, for /3 G M and a sequence of permutations n n ~ Mp >n , the empirical 
permutation measure -+ Mp. This together with Theorem 3.1 can be used to compute the 
limiting density of the degree proportion in a Mallows random permutation. 

7.1. Proof of Theorem 3.8. The proof of Theorem 3.8 has two parts: to show the existence of 
limit of the degree process d n (-) by verifying (6.11) in Corollary 6.3, and the explicit computation 
of the density of the limiting distribution using Theorem 7.1. 

7.1.1. Existence of the Limit. In this section we show that the degree process of a Mallows random 
permutation converges. In light of Proposition 6.1 it suffices to verify that the Mallows random 
permutation is equicontinuous: 
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Lemma 7.1. Let 6 I and ir n ~ Mp, n be a sequence Mallows random permutations. Then 
(7r n )„>i is equicontinuous. 


Proof. Fix e > 0 and 1 < i\ < i^-'-is < n. Let x = ( X\,X 2 , ,x s ) E (0, l] s and y = 

(yi, y 2 ) • • • , y s ) £ (0, l] s and consider the bijection (6.10) between Lt(ii,i 2 ,... ,i s ,x) and fl(ii, * 2 ,..., i s 
If 7r n denotes the image of 7r n under this bijection, then 

1 1 ' 2s 

-|A(7T n ) - A(7f n )| < - V] I \nxa \ - \ny a ] \ < sd + —, 
n n n 

a= 1 

which goes to zero after taking limits as n —> oo and 5 —> 0. Equicontinuity of (7r n ) n >i now follows 
from Corollary 6.3. □ 


The above result and Proposition 6.1 implies that the permutation process 7r n (-) => Wp(-), 
such that for every t > 0, Wp(t) is independent and distributed according to conditional law of 
Qi\Q2 = t, where (Qi, Q 2) ~ Mg. Since the distribution of Mg of (Q 1 , Q2) has uniform marginals, 

rw 

F(Wp(t) <w) = / mp(x,t)dx, (7.3) 

Jo 

and Wf-j(t) has density mp(-,t). Theorem 3.1 then implies that d n {-) => Dp(-), where 


Dp{t) = t + Wp(t) — 2 Mp(t, Wp(t)), 


and Dp{t ) is independent for all t > 0. Therefore, for indices 0 < r\ < r 2 < ■ ■ ■ < r s < 1, 


/ dn{\nr{\) d n {\nr 2 ~\) 
\ n ' n 


^n(rn,r s ]) \ 

n J 


(Di t p, D 2 ) p,... D S 'p), 


where D a> p = Dp(a), for a E [0,1]. 


7.1.2. Calculating the Limiting Density. Fix (3 E M and a E [0,1] and suppose W ~ Wp(a ) be 
distributed as in (7.3). To find the density of D a ,p for a E [0,1], we have to find the density of the 
random variable 


JaA W ) ■= a + W-2Mp(a,W) 

2 /2exp (— \(3{a + 

= a + W + - log I - 


W 


-!)) (sinh(^)sinh(^)) 


- 1 


where 

tpp(a) := e 2 ^ _/3a )csch sinh • 

We begin by establishing properties of the function J a ^p : M —> M defined as J a ,p ( w ) = a + w — 
2 Mp(a, w), for a E [0,1]. Recall that an interval [a, b] is always interpreted as [a V b, a A b] 


Lemma 7.2. Let > 0, a E [0,1], and a c ((3 ) be as defined in Theorem 3.8. Then for J a ^p as 
defined above, the following hold: 

(a) The function J a ^p is strictly convex in M for (3 > 0 and linear for (3 = 0. 






26 


BHASWAR B. BHATTACHARYA AND SUMIT MUKHERJEE 


(b) For a G [0,a c (/3)] the function J a ^ is strictly increasing and for a G [1 — a c (/3), 1], the 
function J a $ is strictly decreasing in [0,1]. 

(c) For a G (a c (/3), 1 — a c (/3)), f/ie function J cij p has a minimum at zq G (0,1), and is strictly 
decreasing in [0, 20 ) and strictly increasing in (zq, 1]. 


Proof. The derivatives of the function J a ^ are 


d 


^a,pi z ) — 7 Ja,p{ z ) — 1 + 


2 <pp(a)e 


-Pz 


and 

J a,/3(-) = 2l PP (°) 


dz u ’ Hy ' ' ipp(a) (1 —e P z ) — 1 

-/3e~ 0z (ipp{a) (l - e~ pz ) - l) - /3e~ 2 P z ipp{a) _ fie^ (tpp{a) - 1) 


(^(a)((l - e P*)- l)) 2 


<Pp(a) ((1 - e P*)- 1) 


2 • 


Note that (p a (P) — 1 = ^jf-i & < 0, for all /3 G M. Therefore, for f3 > 0, J”p(z) > 0 for all zGl, 
and so, J a ^ is a convex function. 

The convexity of J a ^ implies that J' a p is increasing, and J' a ^(z) = 0 has at most one solution 
z 0 in [0,1]: 

z o ■■= jp log ( 1 G (0, !) <=> « e M0), 1 - a c (/3)\, 

where a c (/3) is dehned in Theorem 3.8. Therefore, for a ^ [a c (/3),l — a c (/3)], the function J a ^ is 
strictly monotone: for a G [0, a c (/3)] the function J a p is strictly increasing, and for a G [1 — a c (/3), 1], 
the function J a>( g is strictly decreasing in [0,1]. 

For a G [a c (/3), 1 — a c (/3)], the function J a< p has a minimum at zo, and is strictly decreasing in 
[0, 20 ) and strictly increasing in ( 2 : 0 , !]• □ 


Recall the definition of R(a,/3 ) from (3.10). For notational convenience, for /3 > 0 and a G [0,1] 
define the function h a ,p : [0,1] -A M, 


K,p{ x ) = 


(1 - e~P)y/eP( a+x ) - ehR(a, (5) 

p e \P{a~ x +2) 


._ ■ (7.4) 

y/4eP(eP a - l)(e^ a - eh) + e^ a+x \l - eh) 2 

The above lemma will now be used to complete the proof of Theorem 3.8. Assume that (5 > 0, 
and consider the two cases (recall the definition of J a ,p from (7.4)): 


1: Suppose a ^ [o c (/3),l — a c (/3)]. In this case, J a ,p is increasing in [0,1] (Lemma 7.2) and 
the equation J a ,p(z) = w has a unique solution J~o(w) G [0,1]. Then by the Jacobian 
transformation and direct calculations, the density of J a ,p(W) simplifies to 


9a,p A) 


-j—J n l(w) 
dw a,/3 


= K,A W )- 


The support of J a ,p(W) is [ J a,/?(0), J a ,p( 1)] = [a, 1 - a]. 

2: Suppose a G [a c (/3), 1 — a c (/3)\. In this case, for w G [0,1] the equation J a ,p(z) = w has at 
most two solutions in [0,1] depending on the value of w (Lemma 7.2). 
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2.1: w E [J Q)J g(0), J a ,^(l)] = [a, 1 —a]. This situation is same as the previous case, that is, 
JaA z ) = w ^ ias a llll ' c l lie solution J~p(w) E [0,1], and the density of J a ,/ 3 {W ) simplifies 
to 


9a A w ) 


_d 

dw 


£; J ajH 


M J o}( W )’ a ) = h a,piw)- 


for w E [a, 1 — a]. 

2 . 2 : w £ [J a ,p(0), J a ,p(X)\ = [a, 1 — a]. In this case, J a ,p(z) = w has two solutions given 
by J a,l, i( w ) = -§log7i(w) and Ja,lA w ) = -1 log 72 (w), where 71 ( 10 ) and 72 ( 10 ) are 
roots of the quadratic 


^ + 1 = <pp(a) (l - 7 2 ) =>• cp l 3(a)'y 2 + ye 2 ^ °) + 1 - ^(a) = 0. 


The above quadratic equation is obtained by simplifying the equation J a ,p(z) = w and 
substituting 7 = e~d z l 2 . Then by the Jacobian transformation, the density of J a p(Z) 
is 


9a,p( w ) 



S=1 


(7.5) 


Substituting J~^ /w), J~\ 2 i w ): anf l the density fp(-,a) (7.1), and simplifying (7.5) 
gives g ai p(w) = 2 h a ^{w). Since the function J a ^ has a minimum at zq, and is strictly 
decreasing in [0, zo) and strictly increasing in (zq, 1], the support of J a ,p(W) is 


[Ja,fj{zo), Ja,/3{ 0 ) V J a ,p{ 1 )] — 


1 — a + — log R(a, /3), a V 1 — a 


For the (3 < 0 the result can be proved similarly. Theorem 3.8 holds verbatim even for /5 < 0, if 
every interval [a, b] is interpreted as [a V b, a A b\. 


8. Asymptotics For The Minimum Degree: Proof of Theorem 3.7 

This section gives the proof of the limiting Rayleigh distribution of the minimum degree in a 
uniformly random permutation graph. 

For i E [n] define 

c (i) = l * + 7r "(*) forl<i<^d 

n | 2 (n + 1) — i — 7 T n (i) for < % < n . 

The following lemma shows that the degrees d n (i) can be small (order y/n) only if c n (i) is small, 
which can happen only if i is such that either i or n + 1 — i is small (order y/n). 


Lemma 8.1. For any 7 E (0,oo) 


n 

lim lim P (d n (i) < 7 y/n, c n (i ) > My/n) = 0 . 

M-> 00 n —>00 Z —< 
i=1 

Proof. By symmetry it suffices to show that 

lim lim P (d n (i) < 7 y/n,i + ir n (i) > My/n) = 0, 

M—>00 n—>00 Z —/ 
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which follows if we can show the following: 

lim P (d n (i) < 7 \/n) = 0 , 

71 .—^00 * 


n+1 < - < n+l 

4 — — 2 


lim lim N 

Af-»oo n->oo ^' 


M V™ s'js' n+l 
2 4 


lim lim ) ) 

M->oon->oo ^ z ' 


i) < 7Vn) = o, 


t(*) < 7 \/^, 7 r n(i) = j) = o. 


( 8 . 2 ) 

(8.3) 

(8.4) 


1 M+ri <■,<,. 


Recall by Proposition 5.1, E(d n (i)| 7 r n (i) = j ) = 


_ ■ \ _ (i-l)(n-j)+(j-l)(n-i) _ j(n-2i+l)+(n.+ l)i-2n 


n— 1 


n—1 


Therefore, for < i < z ^, 

E(d n (z)| 7 r n (i) = j) - 7 \/n > j 


n + l 


,) +i (+{H 7v/ " + 2)2 


n+l 

8 ’ 


2 (n — 1 ) 

for all n large enough. An application of Lemma 5.1 now gives 

P (dn(i) < 7 Vn\n n (i) = j ) < 2e'^s. 

On adding over i and j gives 

E P(dn(*) < 7 -\Am + 7 T n (i) > M\fn) < E P(d n (*) ^ 7V / +) < ne“i 28 , 


n + 1 

4 — — 2 


n±lc-ic n +i 
4 - - 2 


from which ( 8 . 2 ) follows. 

Proceeding to prove (8.3), for < i < +1-1, 


E(d n (i)| 7 r n (i) = j ) - 7 y/n >j 
for all M >47 + 8 . Lemma 5.1 gives 


(+VT)) +i - (7V " + 2)a 


i + j 


E(d n (i) < 7V / n|7T ri (z) = j) < 2 E e (* 2 +-? 2 )/ 8 , 

i>i 

which on summing over i and j gives 


E 


/AAJ /’lAJ 

t (i) < 7 \/n) < 2 / e~ x2 ^dx / e~ y2 ^ 8 d 


My/n n+l 

2 - - 4 


’ M /2 


y • 


Since the RHS of the above equation goes to 0 on letting M —> 00 , (8.3) follows. 
Finally, to show (8.4), for 1 < * < and < j < n, note that 

i + j 


E(d n (i)\Tr n (i) = j ) - 7 \/n > j/2 - y^/n - 2 + i > 

for M > 4y + 8. Thus by a similar argument as before, we have 

rM /2 

lim V E F ( d n(*) < 7\/n,7r n («) = j) <2 / 

n—>■ 00 zz ^ / n 

which goes to 0 as M —>• 00 as before. This completes the proof of the lemma 


2 / 32 dz 


a -y 


M /2 


! / 32 dy , 


□ 
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The following lemma now strengthens the above result to show that d n (i) and c n (i) are close for 
those indices i where either i or n + 1 — i is small. 


Lemma 8.2. For any 7 > e > 0 


(a) lim^oo YJi=x F ( d n {i ) < jy/n, c n (i ) > (7 + e)y/n) = 0, 

(b) lim^oc Yli=x > Jy/n, c n (i ) < (7 - e)y/n) = 0. 


Proof. We claim that for every fixed M < 00 , e > 0, 

lim F(\d n (i) - Cn(i)\> Sy/n,Cn(i) < My/n) = 0 . 

n .— 


(8.5) 


1 <i<n 


Proceeding to complete the proof of the lemma using (8.5), note that 
P(dn(*) < 7 y/n, Cn(i) > (7 + e)y/n) 

<P (d n (i) < 7 y/n,c n (i) > Myfn ) +P(| d n (i) - c n (i)\ > £y/n,c n (i) < Myfn). 

Summing over i and letting n —> 00 followed by M —> 00 , the second term goes to 0 by (8.5), and 
the first term goes to 0 by Lemma 8.1. This completes the proof of part (a). The proof of part (b) 
follows by similar calculations. 

Turning to the proof of (8.5), note that by symmetry it suffices to show that 

lim ) 

n—y 00 * 


( i ) - c n (i )| > £y/n, c n (i) < Myfn ) = 0, 


- ra +1 


which follows if we can show 
1 

n— >00 n 


lim- Y' V P(|d n (i) - i - j\ > £y/n\-K n {i) = j) = 0 

i— yoo 71 < < -J 

1 <i<My/n 

To this end, for 1 < i,j < Myfn, 

(*- 1 ){n~j) + {j~ 1 ){n-i) 


\E(d n (i)\TT n (i) =j)-i- j | = 


n — 1 


1 — 3 


< 2M, 


and 


=i) = ( ^ 1) o- 1 >("-7"-^- ) < 2M . 


(n — l) 2 (n — 2) 


4M 2 


Therefore, by Chebyshev’s inequality 

P(|d n (i) - * - j\ > £yfn\ir n {i) = j) < 

This readily gives 

~ Y Y p (\dn(i) ~ i ~ j\ > £y/n\ir n (i) = j) < 


fyfd - 2M) 2 ‘ 


4 M 4 


which goes to 0 on letting n —> 00 , for every M < 00 and e > 0. 


(e^/ra - 2M ) 2 ’ 


□ 
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8.0.3. Completing the proof of Lemma 3. 7. Using the above lemmas we can now complete the proof 
of the theorem. To this end, it suffices to show that for any 7 > 0 

lim P (d n {i) > Jy/n, for all 1 < i < n) = e~ l2 ^ 2 . ( 8 . 6 ) 

n—>■ 00 

Note that 

n 

P (d n (i) > 7 yfn, 1 < i < n) < ^ P (d n {i) > 7 yfn, c n (i) < (7 - e)yfn) + P(c n (?) > (7 - s)yfn, 1 < i < 

1=1 

and so by Lemma 8.2 

limsupP(d n (i) > 7 yfn, 1 < i < n) < limsupP(c n (i) > (7 — e)yfn, 1 < i < n). 

n—> 00 n—> 00 

A similar argument gives 

limsupP(c n (i) > 7 yfn, 1 < i < n) limsupP(d n ('i) > (7 — e)yfn, 1 < i < n), 
n—^oo n—>00 


and so to prove ( 8 . 6 ) and hence the theorem, it suffices to prove the following lemma: 


Lemma 8.3. Let c n (-) be as defined in (8.1). Then for any 7 > 0, 


lim 

n —>00 V 1 <i<n 


( min c n (i) > 


n I = e 7 


(8.7) 


Proof. Note that 


Moreover, 


min Cn(i) > 7 yfn =P (ir n (j) > 7 yfn - j, 1 < j < '/yfn) 
l<.i<.^y/n J 

(n — fyyfn] + 1 ) 

n(n — 1) • • • (n — [_7?iJ + 1) 

(n — I'yy/n] + 1)^”! [n — |_ 7 nJ)! 
n! 


( 8 . 8 ) 


mm c n 
n+1— r Yy/n<j<n 


(j) > 7\/™ 7Tn(*)) 1 < i < 


=P (7Tn(j) < 2(n + 1) - j - 7\/n, n + 1 - 7 ^ < j < n|7T n (i), 1 < i < jy/n) 
(n — 1 — 2 |_7\/ r bl) (n — 2 [7 v^J) • 

( n — L7\A j J )• 


(8.9) 


where the lower bound uses the following argument: The probability of the event is minimized 
when all the Tr n (i), for i e [|_7\AlIL are a t m ost n + 1 — 7 yfn. This minimizes the choices of ir n (j), 
for n + 1 — 7 yfn < j < n. In this case, each i\ n (j ) has (n — 1 — 2 ( 7 \/raJ) choices, and the bound 
follows. 
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Combining (8.9) and ( 8 . 8 ) and taking limits as n —> oo gives the lower bound 
P min c n (i) > 7 y/n 

\l<i<n 

=E I P min c n (j) > 7 y/n vr n (i), 1 < i < 7 y/n ) 1 < min c n (i) > 7 y/n > ) 

\ \n+l—7 x /n<j<n J (l<*<7\/n ) J 

(n — Ltv^J + 1 )L 7 n J(n — 1 — 2 |_ 7 \/raJ) (n — |_7 n J)K n — 2|_7\/™])! 


> 


n!(n — |_7\/nJ)! 


=>- 7 * 


( 8 . 10 ) 


For the upper bound, setting N n := |{1 < i < 7 y/n : n n (i) > n + 1 — 7 y / ro}| and fixing a large 
integer M we have 


( min c n (i ) > 7\/n) <P( min c n (i) > 7 V^, iV n < M) + P(iV n > M) 

l<z<n l<i<n 


<P( min Cn(i) > 7 V™, N n < M) + IE ^, 

l<i<n 1V1 


( 8 . 11 ) 


by Markov’s inequality. Now, since 


\nVn\ 


E N n = Y, 


hVn\ < 2 

— I J 


2=1 


n 


the second term in the RHS of (8.11) to 0 after taking limits as n —> 00 followed by M 7 00 . 
Again, by a similar argument as the lower bound, on the set {N n < M} we have 

P ( min Cn(j) > 7 Vn\ir n {i), 1 < * < 7 V™ ) 

\n+l— 7 y / n< < 7 <n / 

=P (7r re (j) < 2 (n + 1) - j - 7 \/n, ra + 1 - 7 ^ < j < n\ir n (i), 1 <i < 7 y/n) 

(n - 1-2 [ 7 + M) LtAJ ( n _ 2 [ 7 VnJ + M) ! 


< 


(n- /ly/n\)\ 
Therefore, using ( 8 . 8 ), (8.11) and (8.12), 


( 8 . 12 ) 


[ min c n (i) > 7 y/n, N n < M ) 

\l<i<n J 

min c n (j) > 7v / n|7r n (i), 1 <i< 7 y/n ) 1{ min c n (i) > 7 V”}! {An < M} ) 
- 1 —■ r Yy/n<j<n J \<.i<syy/n J 

(n-l- 2 L 7 ^J+M)L^J( n _ 2 L 7 ^J+M)! m ^ . ^ 

- - -:--P mm c n [i) > 7 y/n 

(n-[7V n J) ! Vi<*<7Vh / 


=E P 


< 


—>e 


-7 


(8.13) 


by taking limits as n -A 00 and M — >• 00. This completes the proof of the upper bound, which 
combined with the lower bound ( 8 . 10 ) gives the result. □ 
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